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Abstract 

The coupling of Jupiter's magnetosphere and ionosphere plays a vital role 
in creating its auroral emissions. The strength of these emissions is depen- 
dent on the difference in speed of the rotational flows within Jupiter's high- 
latitude thermosphere and the planet's magnetodisc. Using an azimuthally 
symmetric global circulation model, we have simulated how upstream solar 
wind conditions affect the energy and direction of atmospheric flows. In order 
to simulate the effect of a varying dynamic pressure in the upstream solar 
wind, we calculated three magnetic field profiles representing compressed, 
averaged and expanded 'middle' magnetospheres. These profiles were then 
used to solve for the angular velocity of plasma in the magnetosphere. This 
angular velocity determines the strength of currents flowing between the 
ionosphere and magnetosphere. We examine the influence of variability in 
this current system upon the global winds and energy inputs within the Jo- 
vian thermosphere. We find that the power dissipated by Joule heating and 
ion drag increases by ~190 % and ~185 % from our compressed to expanded 
model respectively. We investigated the effect of exterior boundary condi- 
tions on our models and found that by reducing the radial current at the 
outer edge of the magnetodisc, we also limit the thermosphere 's ability to 
transmit angular momentum to this region. 
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Amongst the eight planets in the solar system, Jupiter, in addition to be- 
ing the largest planet, also has the largest magnetic moment and the largest 
magnetosphere. The magnetosphere interacts with both the solar wind and 
the conducting layer or ionosphere, in the planet's upper atmosphere. These 
interactions can be quite complex and we may use models with some simpli- 
fying assumptions (e.g. axial symmetry) to gain insight into the dynamics of 
the magnetosphere, upper atmosphere and their physical interactions with 
the solar wind. 

Several models of Jupit er's magnetosphere a n d ionosphere have been de- 



veloped in recent studies (iNichols and Cowleyl . l2004l : I Cowley et al.l . l2005l . 



20071 : ISmith and Aylwardl . l2009l ). The se models range frorn deta iled stud- 



ies of the middle magnetosphere onl v (INichols and Cow 



studies of the entire magnetosphere (j Cowley et al 



2005 



2004 ) to global 



2007h and invest! 



gations of the coupled magn etosphere, thermosphere and ionosphere systems 
flSmith and Avlwardl . liooj ) (henceforth SA09). 



One of the important observations for guiding models is the dominance of 
Jupiter's magnetosphere by the rapid planetary rotation. Angular momen- 
tum is transferred from the planet to the disc-like middle magnetosphere via 
ion-neutral collisions in the ionosphere. The magnetospheric plasma exhibits 
a wide range of angular velocities, corresponding to a modest departure from 
rigid corotation with the planet at distances near lo (6-10 Rj) out to regions 
beyond 20 Rj which rotate at ~50 % of the planetarv rate (IMcNutt et al. 



19791 : Imil . I1979I : iHill et all . Il983al : IPontiusl . Il997l : IVasvliunasl . Il983[ ). This 
angular momentum and energy transfer between the ionosphere and magne- 
tosphere is conveyed by two principal current systems. The first of these is 
related to the rotation of the middle magnetosphere. The second is related to 
the interaction of the solar wind with the magnetosp here at the high- latitude 



magnetopause ( IHill et al.l . Il983bl : llsbell et al.l . Il984j ) 



The principal source of plasma for t he middle magnetosphere (' ^20R.t 
to several tens of Rj) is the satellite lo (IBagenal and Sullivanl . Il98ll ) which 
ejects about 500-1000 kg s~^ of sulphur dioxide gas which is then ionised 
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(jKivelson et al.l . |2004J ). logenic plasma initially near corotation will lag fur- 
ther behind corotation as it diffuses radially outwards from the lo torus, 
due to the finitely conducting ionosphere being unable to supply all of the 
necessary angular momentum via the coupling currents. The electric field 
in the neutral atmosphere's rest frame depends on the difference in angu- 
lar velocity between the polar thermosphere and the magnetically conju- 
gate plasma disc, and drives a flow of equatorially directed Pedersen cur- 
rents. Due to current continuity, field-aligned current (FAC) in the steady- 
state must flow both upwards and downwards along the magnetic field lines 
which connect the ionosphere and magnetospheric plasma disc. Downward 
FACs flow from the outermost magnetosphere to the ionosphere. The up- 



the magnetosphere ( 


Cowlev and Bunce. 


2001: 


Hill. 


2001 




Khurana . 


2001; 


Southwood and Kivelson. 


2001 


). These electrons excite 


emissions in the 



upper atmosphere and produce the main auroral oval at ~15 ° co-latitude 



(Satoh et al.. 1996: Clarke et al. 



19991 : iPallier and Prangd 12001 



19981. l2004J:lPrange et al.l . 119981 : IVasavada et al 



Grodent et al 



20031) . Currents flow radially 
outward in the equatorial plane of the magnetosphere and, via the JxS 
force, accelerate the plasma towards corotation. The Pedersen, radial and 
FACs thus represent a complete current 'circuit' coupling the magnetosphere 
and ionosphere. 

The thermospheric angular velocity at Jupiter partly controls the iono- 
spheric Pedersen currents and thus the dynamics of the magnetosphere. We 
do not, howeve r, have many measurem ents of these thermospheric velocities. 
Studies such as lHuang and Hilll (119891 ) and lPontiusI (Il995l ) have attempted to 
model these velocities by coupling the magnetosphere, ionosphere and ther- 
mosphere, with the assumption that angular momentum was transported 
through the thermosphere solely by vertical viscous transport. These stud- 
ies yielded two main conclusions: (i) the relationship between thermospheric 
and magnetospheric angular velocities was a linear one and (ii) thermospheric 
dynamics could be parameterised using an 'effective' ionospheric conductiv- 
ity, which represented the effects of the difference in angular velocity of the 
thermosphere and deep atmosphere (i.e. planetary value). 



Smith and AylwardI (120081 ) and SA09 showed that, for both Saturn and 
Jupiter, meridional advection rather than vertical viscous transfer is the main 
process by which angular momentum is distributed to the high latitude ther- 
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74 mosphere. They also showed that the neutral atmosphere super-corotates, 

75 ultimately due to the sub-corotation of the middle and outer magnetosphere, 

76 at latitudes just equatorwards of the boundary between the field lines of the 

77 middle magnetosphere and the Dungey-Vasyliunas (D-V) layer (region II in 
7B section I3.1.ip . Super-corotation occurs in a region where ion drag forces 

79 (promote sub-corotation) are insignificant compared to Coriolis forces (pro- 

80 mote corotation). These regions, where ion drag tends to zero, lead to the 

81 upwelling of gas which expands and cools adiabatically. This then causes a 

82 pressure gradient which drives poleward flows at altitudes less than 600 km 

83 (pressures higher than 0.04 /xbar). As ion drag is insignificant in this region, 

84 the Coriolis force has no obstruction and can accelerate the gas to super- 

85 corotate. 'Hotspots' were created in these models by converging meridional 

86 winds at the poles while lower latitude regions were cooled. For more de- 
tailed conclusions the reader is referred to SA09. 



87 



89 The modelling of SA09 combin ed the advanced middle magnetosphere 



100 



102 



90 model of iNichols and Cowleyl (120041) with the axisym metric model of the en 



91 tire magnetosphere presented in lCowley et al.l (120051 ). These magnetospheric 



92 inputs were then coupled to a global two-dimensional circulation model of 

93 the Jovian thermosphere. The auroral region in this coupled model is repre- 

94 se nted by the one- dimensiona l auroral thermosphere and ionosphere model 



95 by iGrodent and GerardI (120011 ) . This auroral profile is linearly scaled at each 

96 time step according to the global pattern of auroral conductance (see SA09). 

97 It is this coupled model of the magnetosphere, ionosphere and thermosphere 

98 that we use in the present study of the effects of solar wind variability on 

99 Jupiter's thermospheric flows. There are some necessary minor differences 
between the original SA09 model and the version used herein. These will be 



101 discussed in sections flM and fI77[ 



103 The aim of this study is to see how magnetospheric compressions and 

104 expansions due to changes in solar wind pressure affect steady-state thermo- 

105 spheric flows and temperatures, and the ensuing effect on predicted Jovian 

106 auroral activity. This aspect, as well as the inclusion of a realistic atmospheric 

107 model in our study, implies that this is a natural extension of the previous 

108 studies that we have mentioned. Our basic approach is as follows. We 

109 start with a 'baseline' magnetodisc of equatorial radius, i?MAf=65Rj, where 
no R j=71492 km is taken as the radius of Jupiter. We then produce compressed 
111 and expanded disc configurations (section l2.4p . Using these magnetospheric 
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112 models as input to the atmospheric model, we run for 50 Jovian rotations 

113 until steady-state is reached. 

114 

115 The theoretical background for our study is given in section [2l In sce- 
ne tion [3] we show and discuss our results, in section H] we show cases where our 

117 radial current boundary condition is changed, and we summarise the findings 

118 in section O At present we run our model until steady-state is reached, but 

119 in the future we will aim to simulate the transient, time-dependent effects 

120 caused by relatively rapid variations in so lar wind pressure and thus mag- 

121 netospheric size (e.g. I Cowley et al. (|2003)). These rapid variations in solar 

122 wind pressure cause significant changes in magnetospheric size, and thus 

123 plasma angular velocity, on time scales of 2— 3 hours. On the other hand, 

124 the large inertia of the neutral thermosphere implies that changes in plasma 
126 angular velocity would affect the thermosphere on longer time scales, such 

126 as ~5— 20 hours. This condition allows for the approximation that plasma 

127 angular momentum is conser ved on the shorter t i me sca les associated with 

128 solar wind pressure changes (I Cowley and Buncd . l2003bl ). Time- dependent 

129 effects are usually neglected in most studies of global energy transport in the 

130 Jovian system. Such studies may thus have bearing on the 'energy crisis' at 

131 Jupiter that has evaded a definitive answer for four decades, the fact that 

132 the planet' s exospheric t empe ratures cannot be maintained by solar heating 

133 alone (e.g. iMiller et al.l (120061 )). As a first step, the present study focuses on 

134 the steady state response of the Jovian thermosphere to different magneto- 

135 spheric configurations (i.e. different solar wind pressures). 



136 2. Theoretical Background 



137 
138 
139 
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141 
142 
143 
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In this section we present a summary of some basic theoretical principles 
that we use througho ut thi s stud y. We rely on work tha t has been conducted 
i n pre v ious studies by Hilll (119791) : iPontiud (Il995l . 119971 ) ; iNichols and Cowley 



(120041 ): ICowlev et all f 



20051,120071), and SA09. 



2.1. Ionospheric Currents 

The frictional drag due to ion-neutral collisions within the thermosphere 
causes a lag from corotati on of the thermo s phere that can be represented by 
a 'slippage parameter' K (iHuang and Hilll . Il989l ) 



-KiQj-QM), 



(1) 
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145 or equivalently 



{QT-^M) = ii-K){nj-nM)- (2) 

146 Here Qm, are the angular velocities of the magnetosphere and Jupiter 

147 itself (deep planetary angular velocity) respectively. Qt is the effective rota- 

148 tion angular velocity of the thermosphere (see SA09). K thus represents the 

149 'slippage' of the neutrals from rigid corotation. 



150 
151 



The coupling of the magnetosphere and ionosphere induces an electric 

152 field within the thermosphere's rest frame which then causes ionospheric 

153 currents to flow. The ionospheric height-integrated Pedersen current density, 

154 i p, and the total, a zimuthally integrated form of this current, Ip{9i), are 



155 fICowlev et all (120071 ) . SA09) 

ip =PiT.p {Qt - ^m) Bi, (3) 

166 and 



ipie,)=2np',i:pinT-nM)B,, (4) 

157 where Sp is the height-integrated Pedersen conductance, Bi is the assumed 

158 radial ionospheric magnetic field, 6i is the ionospheric co-latitude, and pi is 

159 the perpendicular distance to the planet's magnetic / rotation axis. 

160 

161 Current continuity requires that there also exists in the magnetodisc a ra- 

162 dial c urrent density, ip, which can also be azimuthally integrated, represented 



163 as Ip (INichols and CowlevI fl2004j ). SA09). We write 



pjp =2piip, (5) 

ip=8nj:pFe{nT-nM), (6) 

164 where pi=Ri sin 9i {Ri is the ionospheric radius), Bi=2Bj {Bj is the equato- 

165 rial magnetic field strength at the planet's surface) and the function Fe{pe) = 

166 Fi{9i) =Bjpf on a magnetic flux shell which intersects the ionosphere at co- 

167 latitude 9i. p^ is the equatorial distance from the planet centre to the field 
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lines lying i n this shell. We a dopt _Bj=426400 nT as Jupiter's dipole equa- 



torial field (IConnerney et al.l . 119981 ) . and -Rj=67350k m as the radius of the 



polar Pedersen layer (radius of the ionosphere) (i Cowley et al.l . 120071 ) . Note 
that the auroral ionosphere is at high latitudes, where the planet's radius is 
~66854km at the 1 bar surface. and Fi are the equatorial and ionospheric 
fiux functions respectively (discussed further in section [231) • The mapping 
between 6i and pe is represented by the equality Fe{pe)=Fi{9i). 

Another result of current continuity with regard to the variation of the 
Pedersen current with latitude is the creation of FACs which fiow from the 
ionosphere to the magnetosphere. The density of these currents (at the 
ionospheric footpoint of the relevant field line) is 



184 
185 
186 
187 
188 
189 
190 
191 
192 
193 



j\\M 



1 dip 
27rR2sin^,- d9, ' 



(7) 



180 where j\\i{Oi) is the FAC density and the sign corresponds to the north- 

181 ern hemisphere whe re the magnetic field points outward from the planet 

182 flCowlev et aLl . l2007h . 



2.2. Effective neutral rotation velocity 

SA09 define Qt as a weighted average of the effective rotation angular 
velocity throughout the thermosphere-ionosphere. In this section we clarify 
what is meant by this, for the sake of completeness. 



Smith and AylwardI (l2008[ l showed that the equatorward current density 
in the ionosphere consists of two contributions: (1) Pedersen current asso- 
ciated with the azimuthal thermospheric velocity and (2) Hall current 
associated with the meridional thermospheric velocity ug. and ue will 
vary, in general, with altitude z in the thermosphere. W e now define a local 
effective angular velocity ojt ( Smith and Aylward ( 2008[ ) . SA09) as follows: 



PiUT = Pi^j + H Ue, (8) 

CTp 

194 where ap and an are the local Pedersen and Hall conductivities respectively. 

195 Integrating over the height of the thermosphere-ionosphere to get the total 

196 equatorward current, we find that Qt can be defined as 
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Spfiy = J apLOrdz, 
197 where Sp is the height-integrated Pedersen conductivity 



(9) 




198 where z is altitude. In these expressions fi^ is a weighted average of the ef- 

199 fective neutral angular velocity ut throughout the thermosphere-ionosphere, 

200 which also contains contributions from meridional winds. 



2.3. Energy Transfer 

In this section we introduce equations which describe the energy transfer 
from planetary rotation to: (i) magne t ospheric ro tation, and (ii) heating of 
the neutral atmosphere. According to Hill ( 200l[ ). the total power per unit 
area of the ionosphere extracted from planetary rotation, P is given by 



P =QjT, (11) 
r =piipBi, (12) 

206 where r is the torque per unit area of the ionosphere exerted by the JxB 

207 force. The smaller component of this total used to accelerate the magneto- 

208 spheric plasma is 



Pm = ^mt. (13) 

209 The remainder of this power is dissipated in the upper atmosphere as heat 

210 and mechanical work 



Pa = (fij - fiM)r. 



(14) 



211 The power Pa consists of two components, as shown by lSmith et al.l ( 120051 ). 

212 One of these is Joule heating, Pj, and the other is ion drag power, P^, which 
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213 is dependent on the sub-corotation of the neutral atmosphere and is then 

214 viscously dissipated as heat. These are given by 



Pj ={Qt - ^M)r, 



(15) 



and 



(16) 



216 
217 



219 
220 
221 
222 
223 
224 
225 
226 
227 
228 
229 
230 
231 
232 
233 
234 



These expressions can then be integrated over the appropriate region of the 
ionosphere to obtain total (global hemispheric) powers. 

2.4- Magnetosphere Model 

The magnetosphere model component use d in this study is e ssentially 
the same as that used by SA09, based on the ICowley et al.l (120051 ) axisym- 
metric model for the entire magne tosphere and the mo r e adv anced middle 
magnetosphere model proposed by iNichols and Cowleyl ( 120041 ). The differ- 
ence between the SA09 model and the one used in this study is that we also 
use the formalism from ICowley et al.l ( 12007! ) to calculate equatorial magnetic 
profiles for compressed and expanded configurations of the magnetosphere. 
Our coupled model requires as input an equatorial profile of magnetic field 
strength, along with the corresponding flux function (the flux function is 
the magnetic flux per radian of azimuth integrated from the given location 
to infinity). For the axisymmetric, poloidal field models which we employ, 
surfaces of constant flux function define a shell of field lines with a common 
equatorial radial distance pe and ionospheric co-latitude 9i. This allows us to 
magnetically map the ionosphere to the equatorial plane using Fi{9i)=Fe{pe) 



(INichols and Cowley 
given by 



20041 ). The ionospheric form of the flux function is 



235 
236 
237 



F, = Bjp^, =BjR^, Sin" 9i. (17) 

The equatorial magnetic field in the middle magnetosphere, Bze, and corre- 
sponding flux function, Fp , in this region are given by the equations below 
(INichols and Cowlevl . Eopj ) 
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238 
239 
240 
241 
242 
243 
244 
245 
246 
247 
248 
249 
250 
251 
252 
253 
254 
255 
256 
257 
258 
259 
260 
261 
262 
263 
264 
265 



FeiPe 



pe 



+ 



2.5peo 



exp 



.2 (P^ 

5' \pec 



Pe_ 

Peo 

5/2' 



5/2' 



. I Rj 
A ( — 

pe 



+ 



m — 2 \ p, 



R.i 



m-2 



(18) 



(19) 



where Bn 



xlO^nT R^, 



:14.501Rj 

-i„-t 



A=5A xlO^ iiT, m=2.71 



oo 



2.841 



3.335 xl0^nT,peo- 

and r(a, z)=JJ^ t"~^e~*(it is the incomplete gamma function. 
These parameters repres e nt an analytical fit to spacecra ft magnetometer data 



fIConnerney et al.l . Il98ll : iKhurana and Kivelsonl . Il993l ). The magnetic field 



model has a grid resolution of 0.01 Rj which, when magnetically mapped 
to the ionosphere, produces footprints of the field lines separated by angles 
equal to or smaller than the thermospheric model's latitudinal grid spacing. 
This is a sufficient condition to sample realistic FAC profiles and thermo- 
spheric fiow patterns within the ionospheric part of the model. 



Using Eqs. (118II19P as a starting point we are able to calculate model 
magnetic fields and fiux functions corresponding to states of differing mag- 
netospheric size. These models are valid within the range of ~5 Rj to near 
the magnetopause, however in this study we employ a rn i ddle m agnetosphere 



with maximum radial distance of 85 Rj. ICowley et al.l (120071 ) assume that 
Jupiter's magnetosphere consists of two components; the middle and the 
outer regions. They take the equatorial magnetic field strength in the outer 
magnetosphere (beyond 65 Rj for their 'baseline' case) to be constant be- 
tween ~5 and ~15nT. Using Eqs. fll8til9p . valid only within the middle 
magnetosphere, we apply their method of compressing and expanding this 
region's magnetic field configuration. We then use our middle magnetosphere 
field model to obtain solutions for plasma angular velocity Qm in this region 
(section 13.1.11) . For the outer magnetosphere we shall use constant, assumed 
values of Qm- 



Us ing the principles of magnetic fiux conservation described by lCowley et al 
(120071 ). we were able to calculate equatorial field profiles for Jupiter's mag- 
netosphere for different values of solar wind dynamic pressure. To compress 
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275 



279 
280 
281 



(resp. expand) a magnetodisc (middle magnetosphere) from an initial radius 
Rmmoj a uniform southward (resp. northward) perturbation field, AB^, is 
applied to our initial m agnetospheric niodel (described by equations Eqs. 



T9|) ). The formalism in lCowley et al.l ( 120071 ) enables us to calculate AB^ as a 



function of magnetodisc radius Rmm- At a given Rmm ('final' disc radius), 
the flux conservation condition is 



-nRlj^AB, = 271 {Fo {Rmm) - Fo (Rmmo)), (20) 

272 where Fq is the initial profile of the flux function (given by Eq. f lT^ ). Rear- 

273 ranging to solve for ABz 



= (21) 



where ABz<0 for a southward field perturbation, and 



AF = Fo (Rmm) — Fq (Rmmo)- (22) 



Using Eqs. fll8tl22p we calculated equatorial magnetic field and flux function 

276 profiles for three different magnetospheric configurations, namely a com- 

277 pressed system, case A with i?MM=45Rj, a baseline system, case B with 

278 -Rmm=65Rj and case C, a n expanded sy s tem w ith i?MM=85Rj. We choose 
-Rmmo=65Rj (as used by Cowley et al. ( 2007f )). These configurations are 



listed in Table [T] and the respective profiles are shown in Fig. [H 



282 Fig. shows how the magnetic field strength varies with equatorial dis- 

283 tance in the magnetodisc for the three cases. The red and green lines show 

284 compressed (case A) and expanded (case C) magnetic field profiles respec- 

285 tively. Case A with a disc radius of -Rmm=45Rj corresponds to a relatively 

286 high solar wind pressure and a strong equatorial magnetic field. Case C, 

287 representing a relatively low solar wind pressure has a magnetodisc radius of 

288 -Rmm=85Rj and a comparatively weak magnetic field. Fig. [Tb shows how 

289 the corresponding flux functions vary with equatorial distance. By definition 

290 the value of the flux function at Pe=RMM has the same value for all cases. 
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2.5. Ionosphere Model 

For simplicity, we use an auroral ionosphere model from the literature 
to derive a global conductivity model. This conductivity model consists of 
both vertical and horizontal variations which we shall briefly summarise in 
sections 12.5.11 and 12.5.21 For a more detailed description of the conductiv- 
ity model employed , the r eade r is referred to the followin K studies, SA09, 
Nichols and CowlevI fl2004 ) and ICrodent and GerardI fl200lh . 



2.5.1. Vertical dependence of conductivity 



The one-dimensional auroral ionosphere model by iGrodent and Gerard 



(I2OOII ) (hereafter GG) developed for Jupiter, is used to establish the altitude 
dependence of ionospheric conductivity. The auroral model uses a two-stream 
electron transport code to calculate auroral electron and ion densities. There 
are two versions of this model, i.e. 'diffuse' and 'discrete' but for our studies, 
both versions produce similar results. Consequently, the diffuse version is 
used, as it covers a greater region of the main auroral oval and polar cap. 

The GG model outputs Pedersen and Hall conductivity profiles for a 
specific thermal structure. However our thermosphere model has a variable 
thermal structure which is a function of latitude. In order to maintain real- 
istic height-integrated conductivities in the model, at each pressure level we 
calculate the conductivity per unit mass as follows (SA09) 



Si = ^, (23) 

312 where i=P or H representing Pedersen or Hall, a is the conductivity and p 

313 is the neutral mass density. Adjacent pressure levels enclose constant masses 

314 of thermospheric gas (hydrostatic equilibrium assumption). Therefore, the 

315 height-integrated Pedersen (Sp) and Hall {J^h) conductivities depend solely 

316 on the profiles of Sj with respect to pressure and not thermal structure. This 

317 leads to a Pedersen conducting layer, which we define as the region with 

318 conductivity greater than 10 % of the Pedersen conductivity at the auroral 

319 ionisation peak, located at pressures of ~0.8— 0.04 phai or at altitudes of 

320 ~350— 600 km above the 1 bar level. 

321 2.5.2. Horizontal conductivity model 

322 The height-integrated conductivity in the inner and middle magneto- 

323 spheres is dependent on the FAG density according to the following equations 
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324 (INichols and Cowleyi . 



20041 ) 



Sp(j||i) = Spo + ^PjUwi 



(24) 



325 where 



326 
327 
328 
329 
330 
331 
332 
333 
334 
335 
336 
337 
338 



340 
341 
342 
343 
344 
345 
346 
347 
348 
349 
350 



0.16j||,+ 

(j||./0.075)2 



2.45 



X 



[1 + (j||i/0.075)2j [1 + expi-Um - 0.22)/0.12)] 



(25) 



where Epo=0. 0275 mho (INichols and Cowleyi . l2004j ) is the background con- 
ductivity due to solar photoionisation, and Epj(j||j) in mho is an auroral en- 
hancement due to the FAC density in /lA m~^. The dep endence of S ^ on 



Kivelson and Russell 



7' l| ,; is c alculated from Eq. using standard formulae (e.g. 
fll995l . p. 201)). The total conductivity in the ionisation region is dominated 
by Ep due to the small values of Sj^. 

In the outer magnetosphere and polar cap regions, conductivity enhance- 
ment is likely to be present sinc e UV a nd IR auroral emissions are detected 
in these regions. ICowley et al.l (120051 ) set Sp=0.2mho (effecti ve Pedersen 



condu ctivity) in these regions in accordance with the theory of llsbell et al. 



(119841 ). To allow for comparison, we employ the same conductivity value in 
these regions. 



2. 6. Coupled Model 

We couple our magnetosphere model with a global numerical model of 
the thermosphere and a global conductivity model of the ionosphere as de- 
scribed in SA09 and section [231 The resolution of the model grid is 0.2° 
in latitude, and 0.4 pressure scale heights in the vertical direction. That 
is, we use pressure as a vertical coordinate, with the lower boundary at 
2 /ibar (300 km above the 1 bar level) and the upper boundary at 0.02 nbar. 
Altitudes are updated in the model assuming hydrostatic equilibrium. For 
simplicity all models are axisymmetric with respect to the planet's axis of 
rotation. This assumption does not greatly influence the basic physics un- 
derlying the conclusions of this study (SA09). It is important to emphasise 
that the assumption of axisymmetry implies zero azimuthal gradients in the 
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351 model variables. This allows us to represent model outputs in two- dimensions 

352 (latitude and altitude) while still using the three-dimensional Navier-Stokes 

353 equations. 



354 



355 Section 6 in SA09 describes the method of coupling the magnetosphere, 

356 thermosphere and ionosphere models. We employ essentially the same 

357 method, with a few minor changes. The same value of the Jovian radius, 

358 -Rj=71492 km is used for both our flux function calculations and atmospheric 

359 modelling. The coupled model in SA09 ran for 200 Jovian rotations to reach 

360 steady-state. Comparisons of height, temperature and azimuthal velocity in 

361 the inertial frame data for case B were made for run-times of 200 and 50 rota- 

362 tions. Calculations show that between both run-times there was a maximum 

363 relative difference of ~0.4 %, ~0.8 % and ~1.2 % for height, temperature and 

364 azimuthal velocity respectively. This difference causes no significant change 

365 in any other parameters obtained from the model and running the model for 

366 50 rotations would save considerable CPU time. Thus for the purposes of 

367 this study, running the model for 50 rotations was considered sufficient to 

368 reach steady state. 

369 2. 7. Solving the coupled equations of thermo spheric and magneto spheric mo- 
no mentum 



Studies such as iHilll (jl979[ ) and iPontiud (119971 ) have shown that for the 



middle magnetosphere to be in a steady state, the radial gradient of the out- 
ward angular momentum flux of iogenic plasma must be equal in magnitude 
to the torque per unit radial distance on that plasma. The plasma model 
that describes the middle magnetosphere is based on four equations 



Pe dpe ^ M 

= -^nh^ P^^^ - ^^)]' (27) 

Pel-Dzel dPe 

Ep = Sp(j||,), (28) 

nT = nT{nM,^p), (29) 

376 where M=1000kg s~^ is the assumed mass outflow rate from the lo torus 

377 and is the upward FAC density in the ionosphere. 
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379 These equations describe the inter-dependence of magnetospheric angular 

380 momentum per unit mass {pIQm), FAC density (j||i) and Pedersen conduc- 
tance (Sp). Eq. fl29l) represents the output from the thermospheric model 
component, which is forced by raagnetospheric inputs of QM^Pe)- Eq- 



381 
382 



394 



396 



400 



410 
411 
412 



383 is the Hill- Pontius equation (IHilll . Il979l : iPontiud . 119971 ) with a modification 



384 by SA09 to include effects of neutral thermosphere flow, represented by Qt- 

385 This equation balances torques caused by the outward diffusion of the disc 

386 plasma and the JxB force associated with the magnetosphere-ionosphere 

387 coupling currents. Eq. (1271) is used to calculate the FAC in the ionosphere. 

388 An increase in field-aligned current should have an effect on angular velocities 

389 ^MiPe), through enhancement of the ionospheric conductivities. We account 

390 for this in Eq. (1251) . representing Eqs. ([M] and [22]) , which describes how en- 

391 hancements in also affect the Pedersen conductance Sp (higher flux of 

392 precipitating auroral electrons increases the production rate of ionospheric 

393 plasma) . 



395 Our method for s olving these equations is the same as that in SA09 and 



Nichols and Cowleyl (120041 ) (who originated this model). This is essentially a 

397 shooting method which varies the value of Qm at the outer edge of the disc 

398 until the solution, integrated inwards from this location, smoothly joins an 

399 appropriate 'inner disc' analytical solution at 12 Rj. We set the azimuthally 
integrated radial current at the outer edg e of the disc to a value of 10 MA as 



401 our outer boundary condition (following iNichols and Cowleyl (120041 )). whilst 

402 we have near-rigid corotation of plasma as an inner boundary condition. 

403 We need, however, to ensure that the height-integrated Pedersen conduc- 

404 tivities at the poleward ionospheric boundary of the magnetodisc field line 

405 region, Sp(disc), and at the equatorward boundary of the outer magneto- 

406 sphere region, Sp (outer), join smoothly together to avoid discontinuities at 

407 this interface. This is particularly important for large compressions such as 

408 that of case A. A Gaussian function was used to extrapolate Ep from the 

409 magnetodisc into the outer magnetosphere region. We ensured the Gaussian 
function would terminate with a polar value equal to the chosen background 
Ep in the outer magnetosphere, and that this transition would occur with 
a small latitudinal scale (0.2°). The amplitude and centre of the Gaussian 

413 function were calculated using the gradient of dTip/dO at the poleward edge 

414 of the disc region. We further discuss the resulting profiles of Ep in sec- 

415 tion [3X2l 
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416 3. Results and Discussion 



418 
419 



422 



424 
425 
426 
427 
428 
429 
430 
431 
432 



434 
435 
436 
437 
438 
439 
440 



In this section we present tlie results obtained from our modelling. We 
firstly discuss results concerning angular velocities, conductivities and cur- 
rents. Then we proceed to discuss the thermospheric fiows and energies. 



420 3.1. Angular velocities, conductivities and currents 

421 3.1.1. Angular velocities 



Combining models from iNichols and Cowleyl (120041 ) and I Cowley et al 



423 ( 120051 ) of the middle and outer magnetospheres, one can essentially divide the 
entire magnetosphere into four regions, labelled I-IV. The Dungey-type inter- 
action of the polar open field lines with the solar wind takes place in region 
I. The closed field lines of the outer magnetosphere are involved in Dungey 
and Vasyliunas cycles (associated with mass loss from the disc) in region II. 
Regions III (shaded region in figures) and IV represent the middle magne- 
tosphere (magnetodisc) and the corotating inner magnetosphere respectively 
(see Fig. 12]). As stated in section [2^ VLt is a weighted average of the effective 
angular velocity throughout the thermosphere-ionosphere, computed over all 
altitudes at each co -latitude 6*,. VLm in region I has a constant value of ~0.1 r2j 



433 ( llsbell et al.l . Il984[ ). Region II also has a fixed value oi VtM that depends o n 



magnetospheric size, in accordance with observations ( ICowley et al.L 120071 ). 
The profiles of VLm in regions I, II and III are joined smoothly across their 
boundaries with the use of hyperbolic tangent functions. The plasma angular 
velocity profiles for regions III and IV are calculated using Eqs. (1261- 129|) by 
the model. 



Fig. |2] shows how the thermospheric (solid lines) and magnetospheric 

441 (dashed lines) angular velocities vary in Jupiter's high latitude region for 

442 our three cases. We also show the region boundaries used in our model 

443 and the magnetically mapped location of lo in the ionosphere. Case B, our 

444 'baseline' is shown in blue. At low latitudes, rigid corotation with Jupiter's 

445 deep atmosphere is maintained. At the higher latitudes (> 60 °) the mag- 

446 netosphere (represented by VLm) sub-corotates to a greater degree than the 

447 thermosphere (expressed by VLt)- The shape of these VLm and VLt profiles are 

448 similar to those obtained in the studies of SA09. VLm and VLt profiles for case 

449 C, our expanded case, are represented by green lines. These profiles resemble 
those of case B but they possess slightly smaller angular velocities in region 



450 
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II. For case A, and ^It are shown by the red lines. Both Q,m and Q,t indi- 
cate sub-corotation to a lesser extent tha n the respective pro files from cases 
B and C, in agreement with the study of ICowley et al.l (120071 ) who modelled 
Qm, assuming simplified profiles for Qt (where K=0.5 ). We thus show that 
the thermosphere and magnetosphere for compressed configurations corotate 
to a greater degree than in the case of expanded configurations. Our plotted 
profiles quantify this result for both Qm and Qt- 

3.1.2. Conductivities and Currents 

Previous studies of the effect of solar wind-induced compressions and ex- 
pansions of Jupiter's magnetosphere have shown that magnetospheric com- 
pressions reduce ionospheric and parallel currents (in the steady state). Ex- 
pansions on the other hand, have the opposite effect due to the increased 
transport of angular mornentum t o the magnetosphere (jSouthwood and Kivelsonl . 
200ll : ICowlev and Buncd . l2003al : ICowlev et all l2007h . Our profiles in Fig. [2] 
confirm and quantify the expected angular velocity profiles of both the ther- 
mosphere and magnetospheric plasma in the steady state, when the rate of 
addition of angular momentum to the plasma (at a given radial distance), 
due to the magnetosphere-ionosphere currents, exactly balances the rate of 
removal due to the radial plasma outfiow. We consider the solutions for VLm 
and VLt in more detail in section I3.1.1[ The weaker average magnetic field for 
the expanded cases, combined with the finite ionospheric conductivity, leads 
to lower VLm values, despite increased rates of angular momentum transport 
in the system. In this section we present our quantitative findings regarding 
ionospheric conductivities and currents for the different magnetospheric con- 
figurations of our coupled system. 



The variation of height-integrated true Pedersen conductivity Sp for our 
three magnetospheric cases is shown in Fig. [3^, where cases A-C are repre- 
sented by red, blue and green lines respectively. The magnetically mapped 
location of lo in the ionosphere is shown by the black dot and the mag- 
netospheric regions used in this study are marked and separated by black 
dotted lines. All three cases have peaks just equatorward of the region III 
/ II boundary — characteristic features of the VLm solutions (Eqs. (12611291) ) 
— and then fall to the assumed conductivity value in regions II and I. Cases 
B and C have similar profiles and peak values close to those calculated in 
SA09, whilst case A has a peak that is significant ly higher than both of these 
cases. The profile for case A resembles that from lNichols and Cowleyl (120041 ) 
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490 



for the near-rigid corotation approximation where (1 — Q^j/Qj) << 1, which 
are conditions met by case A in regions IV and III. Another feature that dis- 
tinguishes case A is that the peak conductivity is shifted poleward shghtly 

491 compared to cases B and C. This is partly due to the model method which 

492 connects the Pedersen conductivity in region III with the fixed value in re- 

493 gion II for case A (see section [22]). The poleward shift is also due to the 

494 higher S p required in case A in order to achieve the prescribed value of radial 

495 current at the outer edge of the magnetodisc (poleward boundary of region 
III) (see section [2T|l . 



496 
497 
498 



501 
502 
503 



Fig. [nb shows how the slippage parameter K varies with latitude for 

499 our three magnetospheric cases. The profiles for K indicate the ratio be- 

500 tween thermospheric and magnetospheric angular velocities with respect to 
Jupiter's planetary rotation velocity {K={flj — Qt) / (^j — ^m))- Positive 
values for K represent situations when both the thermosphere and magne- 
tosphere are sub-corotating or super-corotating with respect to the planet, 

504 as seen in region IV, II and I. Negative K values represent situations where 

505 the thermosphere and magnetosphere are undergoing opposing motions i.e 

506 one is super-corotating whilst the other is sub-corotating. This is seen just 
equatorward of lo's magnetic footprint on the ionosphere and for the latitu- 
dinal majority of region III. This distinction is important because the last 

509 half degree of latitude in region III maps to the largest part of the equatorial 

510 magnetosphere. 



507 
508 



511 



512 Fig. shows the corresponding variation of azimuthally-integrated Ped- 

513 ersen current with latitude. The colour code is the same as that in Fig. |2j 

514 Profiles for cases B an d C follow a s i milar trend to the steady-state Peder- 



518 



515 sen current profiles in ICowley et al.l ( 120071 ) , whilst the profile for case A is 



516 different within regions III and II, due to conditions comparatively nearer to 

517 corotation. 



519 Fig. HJd shows the azimuthally integrated radial currents through the mag- 

520 netospheric equator for all three cases in regions IV and III and how they 

521 vary with radial distance. The radial currents for cases B and C show a 

522 's-curve' structure which is consistent with previous studies such as SA09. 

523 Case A however, shows a more linear relation between the equatorial radial 

524 distance and azimuthally integrated radial current which is not seen in the 

525 more expanded case of SA09 but is consistent with the near-rigid corotation 
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526 approximation conditions presented in iNichols and Cowleyl (120041 ) . As previ 



527 



ously noted, this near-rigid corotation condition applies to case A throughout 

528 regions IV and III. We also note that, as mentioned in section [2771 our outer 

529 boundary condition is that the radial current value at the region III / II 

530 boundary is 100 MA. The case A curve in Fig. |3)d does not quite reach 

531 this value due to the joining of the Pedersen conductivity across regions III 

532 and II (see section . A hyperbolic tangent function is used to smoothly 

533 join the Perdesen conductivity across regions II and III, using information 

534 from a few points either side of this boundary. This leads to a smoothing of 

535 the disc solution near its outer edge, leading to a slightly different value of 

536 the azimuthally integrated radial current at this location. This curve does 

537 demonstrate however, that Ip and Ip in case A have to increase very rapidly 

538 in the outer magnetodisc in order to satisfy the boundary condition. Since 

539 there is no a priori reason why /poo should be independent of magnetosphere 

540 size, we will also investigate, later, the effect of varying the boundary condi- 

541 tion upon the resulting profiles of current and angular velocity (section H]). 

542 

543 FAC densities are plotted against latitude in Fig. O For all three cases, 

544 FAC densities have three positive peaks, the first two lying on either side of 

545 the region III / II boundary and the third lying on the region II / I boundary. 

546 Positive peaks correspond to upward directed FACs that produce aurorae. At 

547 the boundary between region III and II, the negative peaks indicate strong 

548 downward-directed FACs whose magnitude is dependent on the equatorial 

549 radius of region III (Rmm)- The main auroral oval is represented by the 

550 peak at ~73 ° latitude. Our model suggests that there would also be weaker 
651 more distributed aurorae poleward of the main oval, represented by the sec- 

552 ond and third peaks at ~75 ° and ~80 ° respectively. A relatively dark region 

553 would arise from the trough at ~74 ° latitude, cre ating 'dark rings'. The lat- 



554 ter feature is also obtained in previous studies by ICowley et al.l (120051 . 120071 ) 

555 but at present, we lack the observations required to constrain the value of 

556 downward. The strong downward FACs at ~74° are due to the significant 

557 changes in Pedersen current on crossing the boundary between regions II and 

558 III, which in turn is due to the changes in magnetospheric and thermospheric 

559 angular velocities. The Pedersen conductivity in the model also changes sig- 

560 nificantly across this boundary, which also contributes to a large magnitude 

561 for jiu. The strongest downward FACs in our calculations are eve n less con 



562 strained by observations, but they also occur in the modelling of iTao et al 



563 (120091 ) who also used a coupled magnetosphere-thermosphere approach. 



19 



564 



665 Our ca lculations shown in FiRs. [ 2 MI all support the expected trends de- 



566 scribed in ISouthwood and KivelsonI (j200l[ ). The angular velocity profiles 



667 (Fig. [2]) for both the thermosphere and magnetodisc show that there is a 

568 greater degree of sub-corotation for more expanded magnetospheres, corre- 

669 sponding to lower solar wind dynamic pressures. This is due to the ther- 

570 mosphere being able to transfer momentum to a compressed magnetosphere 

571 (stronger field) with greater efficiency than a larger, expanded one. The 

572 Pedersen conductivities (Fig. [3ti), FAC densities (Fig. |5]) and azimuthally- 

573 integrated Pedersen and radial currents (Figs. H^-b) all show an increase in 
674 region III (shaded) for expanded magnetospheres. In this region, the inte- 

575 grated auroral FAC for case A is ~50— 60 % of that cases B and C suggesting 

576 that auroral emission would be greater for an expanded magnetosphere than 

577 a compressed one. Our currents naturally have similar values to those ob- 

578 tained in SAO 9. They also show sirn ilar t rends and profiles to studies such 



579 as those of ICowley et al.l ( l2005l . 120071 ) and iTao et al.l ( 120091 ). Our study is an 



580 extension of these works in the sense that we use an atmospheric circulation 

581 model coupled to three distinct magnetospheric configurations. 

582 3.2. Thermo spheric flows and energies 

583 Fig. E] shows momentum balances for our compressed and expanded con- 

584 figurations in both the low and high altitude regions. Fig. [7] shows the ther- 

585 mospheric fiows, temperature distributions and power dissipated per unit 

586 area for all three model configurations. Results for each case are displayed 

587 in the columns of the figure. 

588 3.2.1. Thermo spheric flows 

589 According to SA09, meridional advection is the main process by which 

590 angular momentum is transferred to the high latitude thermosphere. Ad- 

591 vection (combination of the horizontal and vertical advection of momentum 

592 by winds blowing along and across fixed pressure surfaces) is just one of the 

593 means by which momentum is changed at a fixed location within the ther- 

594 mosphere. In Fig. E]we present force balance diagrams at low (a-b) and high 

595 (c-d) altitudes for cases A and C. The force colour codes are in the figure 

596 caption. Considering the high altitude region first, advection and other zonal 

597 force components (ion drag and Coriolis) are small. Thus, the pressure gra- 

598 dient is balanced almost perfectly by the Coriolis force. This force balance 

599 creates a sub-corotational fiow with a small equatorward component. We 
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600 now consider low altitudes near the Pedersen conductivity peak, where the 



601 



ion drag term Jx B is strong. Coriolis, pressure gradient and ion drag forces 

602 are not balanced. Thus, a significant advection term arises to restore equi- 

603 librium, resulting in a region of strong poleward acceleration (see Figs. [31- f). 

604 The resulting meridional fiow at low altitudes is thus polewards and trans- 

605 ports heat to the polar region. 



606 
607 
608 



Figs. [7^-c show how the thermospheric azimuthal velocity in the coro- 
tating reference frame varies within the high latitude region for the different 
609 cases. Positive (resp. negative) values of neutral azimuthal velocity indi- 
cate super (resp. sub) -corotating regions. Arrows indicate the direction 
of meridional flow, and the white line the locus of rigid corotation. The 
magnetospheric region boundaries are plotted with the dotted black lines. 

613 We can see a broad azimuthal jet (blue area) in regions I and II that sub- 

614 corotates to a greater degree with an increase in magnetospheric size. Also 

615 present is a super-corotational jet (dark red region) just equatorward of the 
region III / II boundary, visible in Fig. |2J Ion drag (see Fig. [H]) gives rise 
to the sub-corotational azimuthal flows seen in regions I, II and III. As the 
magnetosphere expands, the JxB term increases and azimuthal flows sub- 

619 corotate to an even greater degree. Advection forces arise due to the lack 

620 of equilibrium at low altitudes, causing an accelerated poleward flow whose 

621 velocity increases by ~90 % from case A to C. The effect of advection can be 

622 seen in Figs. Eli-f, which show meridional flows in the high latitude region. 

623 This accelerated flow transports energy from Jo ule heating, deposi ting it at 



610 
611 
612 



616 
617 
618 



624 



632 



higher latitudes and forming a polar 'hot spot' ( ISmith et al.l . 120071 ). Super- 

625 corotation occurs at latitudes where zonal ion drag and advection forces are 

626 negligible compared to the Coriolis force, which can then accelerate the flow 

627 beyond corotation. At high altitudes, forces are essentially balanced. Thus 

628 high altitude zonal flows now have an equatorward component. Therefore, 

629 meridional flows show a poleward low-altitude flow a nd an equato r ward high- 



630 altitude flow consistent with the previous studies of ISmith et al.l (120071 ) and 

631 SA09. 



633 Figs. Uk-^ show thermospheric temperature distributions. The tempera- 

634 ture scale is shown on the colour bar. Magenta and solid grey contours enclose 

635 areas where Joule heating and ion drag energy inputs exceed 20 Wkg~^ and 

636 dashed grey contours highlight regions where ion drag decreases kinetic en- 

637 ergy at a rate greater than 20Wkg~^. A uniform rate of 20Wkg~^ gives 
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638 



an integrated energy input rate of the order of lOOmWm"^ (1 ergscm"^ 
639 =lmWm~^) within the Pedersen conducting layer. This integrated rate is 
of a similar order of m agnitude to the estimated total IR auroral emission 



640 
641 



(200 ergs cm ^s ^ from iDrossart et al.l ( ll993l )). We see significant energy in 



642 put from Joule heating and ion drag at low altitudes and between 73 °— 85 ° 

643 latitude due to their dependence on the strength of the current density j. 

644 j is proportional to the difference between local thermospheric and plasma 
646 angular velocity, and to the Pedersen conductivity ap which peaks at low 

646 altitudes. There lies a narrow region of high altitude Joule heating just 

647 equatorward of the region III / II boundary in expanded cases due to the 

648 large shear between and Qm- The decrease in kinetic energy (grey dashed 

649 lines) occurs as the ion drag force now acts to accelerate thermospheric flows 

650 towards corotation (see Figs, [nt-d where ion drag is eastwards). The remain- 

651 ing feature of prominence is the large 'hot spot' at low altitudes in region I 

652 as discussed above. The peak temperature of the 'hot spots' increase from 

653 ~560 K in case A to ~695 K in C. 



654 3.2.2. Atmospheric and magneto spheric energy 'budget' 

655 Figs. [7|-1 present the total, magnetospheric, atmospheric, ion drag and 

656 Joule heating power per unit area (see Eqs. flimi6l) ) for each magnetospheric 

657 configuration. The colour code indicates the form of energy dissipation (blue 

658 curve shows total power). Total power is the sum of the magnetospheric 

659 and atmospheric powers and atmospheric power is the sum of Joule heating 

660 and ion drag powers. Integrating the power per unit area over the model 

661 hemisphere gives us the global powers for each of these mechanisms, shown 

662 in Fig. [HI 

663 

664 Figs. [7|-1 exhibit peaks in power per unit area just equatorward of the 

665 interface between regions III and II due to the large ionospheric current asso- 

666 dated with the breakdown in corotation of the magnetodisc (region III). Most 

667 of the energy dissipated in region III accelerates the magnetospheric plasma 

668 towards corotation. This magnetospheric power dominance diminishes for 

669 the expanded magnetosphere, in which more power is dissipated in the at- 

670 mosphere via Joule heating and ion drag. Region II is dominated by these 

671 atmospheric power terms in cases B and C whilst the same region in case 

672 A is still noticeably magnetospheric power-dominated. Atmospheric power 

673 is the major form of energy dissipation in region I. The atmospheric power 

674 dominance in regions II and I is mainly due to the low assumed values for 
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675 magnetospheric angular velocity Qm (see Fig. [2]). The difference {Qt — ^m) 

676 is largest in these regions, which produces relatively large ionospheric Peder- 

677 sen currents and atmospheric power. The low value of Qm produces a mag- 

678 netospheric power that remains low compared to other regions. Note that 

679 relative amounts of energy provided to the atmosphere and magnetosphere 

680 on any flux shell depends on the difference {Qt — ^m) through Eqs. 
They are equal when flM=0-5ftj. 



681 
682 



687 



706 



708 



710 



683 These figures suggest that as Jupiter's magnetosphere is compressed a 

684 higher proportion of the total power of planetary rotation (in the steady 

685 state) is transferred to the magnetosphere via the magnetic field, and by 

686 contrast, as the magnetosphere is expanded more heat is dissipated in the 
atmosphere. 



689 Fig. [8] shows the integrated ion drag (blue). Joule heating (green) and 

690 magnetospheric (red) power per hemisphere for each case and how these 

691 powers are distributed in the open and closed field line regions. Powers in 

692 the closed field regions lie below the dashed white line whilst powers in the 

693 open field regions lie above it. Integrated powers within region I remain 

694 essentially unchanged for both atmospheric (ion drag plus Joule heating) 

695 and magnetospheric mechanisms due to the assumed constant value of the 

696 magnetospheric angular velocity for all cases. Atmospheric power increases 

697 significantly with magnetospheric size, by a factor of ~3 from case A to 

698 C. Magnetospheric power shows a slight decrease between cases A and C 

699 and is a maximum in case B. Magnetospheric power is proportional to the 

700 torque, which increases with magnetospheric expansion, and the magneto- 

701 spheric angular velocity which decreases with expansion. One would thus 

702 expect that, given a large enough range of magnetospheric sizes, the magne- 

703 tospheric power would display a non-monotonic profile that is 'modulated' 

704 by the difference between the angular velocity of the neutral thermosphere 

705 and plasmadisc. 



707 We now compare our int egrated powers per hemisphere with those calcu- 



lated by lCowley et al.l (120071 ) to take into account how using a two-dimensional 



709 model of Jupiter's thermosphere changes previous theoretical power esti- 
mates. We only compare case s A an d B with the 'intermediate' and 'base- 



711 line' cases from ICowley et al.l (120071 ) respectively, as they have comparable 



712 magnetodisc radii. The total integrated magnetospheric powers per hemi- 
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733 
734 
735 
736 
737 



744 



sphere are compara ble in magnitude; 'intermediate' and 'baseline' cases in 



Cowley et al.l (120071 ) are ~85 % those of our cases A and B. The difference be- 
tween atmospheric powers is greater, since this study uses a detailed thermo- 
spheric model . In re gion I, cases A and B have atmospheric powers ~35 % of 



Cowley et al.l (120071 ) 's 'intermediate' and 'baseline' cas e s. In regions III and 



718 II (closed field) the atmospheric powers in I Cowley et al.l (120071 ) are ~70— 80 % 

719 those in cases A and B. Thus, within the closed field region of cases A and B, 

720 the inclusion of a detailed thermospheric flow model has led to more energy 

721 transferred from the thermosphere to the magnetosphere for accelerating the 

722 magnetospheric plasma and more heat dissipated within the thermosphere 

723 via Joule heating and ion drag. 

724 4. Effect of outer boundary conditions 

725 The results for case A exhibit a relatively large peak for the Pedersen 
conductivity and FAC density just equatorward of the boundary between 



727 region s III and II (Fig. [5]). Previous studies such as lSouthwood and Kivelson 



728 (I2OOII ) suggest that such peaks for compressed magnetospheres should be 



729 smaller in magnitude than those for more expanded magnetospheres. This 

730 is in contrast to what we observe in section [31 If the radial current at the 

731 region III / II boundary, Ip^o is larger than realistic values for a compressed 

732 magnetosphere, large FACs would develop in the poleward part of region 
III to satisfy the Ipoo boundary condition. In this section we present model 
outputs with smaller values of I poo- We select illustrative I poo values for each 
case which minimise the variance in the current profiles throughout regions 
III and II. Decreasing Ipoo decreases current gradients in the well constrained 
Qm model for region III and allows for a smoother transition to region II, 

738 whose Qm profile is poorly constrained (due to paucity of observations). 

739 4-^- Outer boundary conditions for a compressed magnetosphere 

740 To commence this part of our study, we ran case A but with /poo =45 MA, 

741 to see whether any significant changes would arise in the currents at the 

742 region III / II boundary. These results are shown in Fig. [9l where blue and 

743 red lines represent cases Aioo and A45 respectively (subscript denotes the 



value of Ipoo] 



745 /^.i.i. Comparison of angular velocities and currents for case A 

746 In this section we compare angular velocities and current-related param- 

747 eters for cases A45 and Aioo (case A in section [3]). The differences between 
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748 these cases essentially lie within 72— 78 ° latitude and so our subsequent dis- 

749 cussions will focus on this range. 



750 



754 
755 



751 Fig. [9^ shows the influence of I poo on the VLm and VLt values. Case Aioo 

752 with the higher /poo value also has a higher torque on the disc plasma, which 

753 results in a smaller difference VLt — VLm near the disc boundary (74—75°). 
Equatorward of 74°, the VLt and VLm profiles are very similar for both cases. 
The outer disc region 74—75° thus develops a strong FAC in case Aioo in 

756 order to satisfy the higher Ip^o imposed. 

757 

758 Due to the smaller values of VLm and radial current for case A45, Pedersen 

759 conductivity values (Fig. IHb) are significantly smaller near the region III / II 
boundary compared to case Aioo- The different gradients in VLm and VLt for 
case A45 cause a slight equatorward shift in the conductivity peak compared 
to Aioo- The 'slippage' parameter for A45 in Fig. [Qt differs only slightly from 
that of case Aioo due to its smaller VLt — differences. 
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The azimuthally-integrated radial, Pedersen and FACs for cases A45 (red 
line) and Aioo (blue line) are shown in Figs. [9ll-f respectively. The other labels 
767 are the same as in Fig. HI For Aioo the magnetos phere near-rigidly corotate s 



768 
769 
770 



774 
775 
776 
777 
778 
779 

780 



with ^Ij throughout regions IV and most of III (INichols and Cowleyi . |2004J ). 
The A45 radial current profile resembles those for expanded cases, due to 
the magnetosphere sub-corotating to a greater degree (see Fig. The 

771 Pedersen current for case A45 has a smooth, almost linear transition across 

772 and through regions III and II as opposed to the abrupt cutoff at the region 

773 III / II boundary in Aioo. The A45 profil e quantitative l y rese mbles Pedersen 



currents for expanded cases and those in I Cowley et al.l (120071 ). For A45, FAC 
profiles are similar to those for Aioo with the exception that the magnitude 
of the peaks just equatorward of the region III / II boundary and the trough 
are significantly smaller. For lower Jpoo values, then, the main auroral oval 
would be significantly dimmer, and possibly more similar to the putatively 
weak auroral signature at the region I / II boundary. 



4- 1.2. Thermo spheric distributions for case A 

781 Here we discuss the changes in azimuthal and meridional velocity as well 

782 as the temperature distribution, which arise from setting Jpoo=45 MA for case 

783 A. All conventions and colours in Figs. [9^-h are the same as those in Fig. [71 



784 
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785 Figs, ink-h show the distribution of azimuthal and meridional velocities 

786 across the high latitude region for case A45. In Aioo (Fig. Uh) there is a 

787 large sub-corotational jet in regions II and I with the strongest degree of 

788 sub-corotation just poleward of the region II / I boundary. In A45, the large 

789 sub-corotating jet now has two regions of strong sub-corotation, the new one 

790 being just poleward of the region III / II boundary. These two strong sub- 

791 corotational jets within the larger jet are evident in Fig. [H^ at the region 

792 boundaries (dotted black lines), where there are large changes in magneto- 

793 spheric angular velocity. Meridional velocities for case A45 follow the same 

794 trend as in Aioo (Fig. [Hi) where there is a poleward flow at low altitudes 

795 and an opposite flow at high altitudes. The main difference between the 

796 meridional flows is that localised accelerated flows (high altitude in region 

797 III and low altitude in regions II and I) have larger velocities in case A45 due 

798 to larger advection terms. 
799 

800 Fig. [9} shows the temperature distribution for case A45. Comparing this 

801 with case Aioo (Fig- Ek) indicates that energy input via Joule heating (ma- 

802 genta contours) and ion drag (solid grey contours) is greater in A45. The 

803 larger energy input, predominantly in region II, is caused by larger shear 

804 between thermospheric and magnetospheric angular velocities. This leads 

805 to a slight increase in thermospheric temperature (~6 %) most evident in 

806 region I, the polar 'hotspot' into which auroral heat energy is transported by 

807 meridional winds. 

808 4-^- Outer boundary conditions for the baseline magneto sphere 

809 For our baseline case, case B, the smallest variance in current profiles 

810 occurred with /poo=68MA. We compare this case Bgg with the original Bioo 

811 case in Fig. dDi 

812 i. Comparison of angular velocities and currents for case B 

813 As for the compressed magnetosphere, we compare angular velocities and 

814 current-related parameters of the Bqs and Bioo models in the 72—75 ° latitude 

815 range where significant differences arise. Fig. [TUb compares the variation of 

816 magnetospheric and thermospheric angular velocities for cases Bgs and Bioo. 

817 In region III, both the magnetosphere and thermosphere for Bqs are slightly 

818 sub-corotating compared to the Bioo- The Pedersen conductivity and 'slip- 

819 page' parameter for Bgs also have similar profiles, but with smaller mag- 

820 nitudes in region III, to those for Bioo. These minor differences are caused 
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821 by smaller thermospheric and magnetospheric angular velocities in region III. 

822 

823 Fig. [TOld-e show azimuthally-integrated radial and Pedersen currents. The 

824 corresponding FAC density as a function of latitude is shown in Fig. [TOl . 
826 The radial current profile is smaller in magnitude than that of Bioo- The 

826 Bgs Pedersen current has a single small peak at the region III / II boundary 

827 in contrast with the sharp peak in Bioo- The FAC density has a smaller, 

828 slightly broader peak in region III, suggesting a low intensity auroral oval 

829 compared to the Bioo case. The absence of strong downward FAC for Bgs 

830 also suggests that the method used to join the region III currents with the 

831 region II currents could produce artefacts for relatively large values of Jpoo- 

832 

833 As with case A, the fine structure around the boundary between the mid- 
834 die and outer magnetosphere in Figs. [TOb-f has been removed by decreasing 

835 the value of /poo from 100 MA to 68 MA for our baseline case. 

836 J^.2.2. Thermospheric distributions for case B 

837 Here we discuss the minor changes in Jupiter's thermospheric flows and 

838 temperature distribution made by changing the radial current boundary con- 

839 dition value for our baseline case. All conventions and colours are the same 

840 as those in Fig. [71 

841 

842 The azimuthal velocity in the high latitude region is shown in Fig. [TI% . 

843 We expect a slight increase in sub-corotation throughout region III (see 

844 Fig. [TOk) compared to Bioo- This is evident by comparing Fig. [TO^ with 

845 Fig. [7^ where we can see that the region of super-corotation (dark red) has 

846 diminished for Bgs- The meridional velocity distribution is shown in Fig.fTOh. 

847 The high altitude localised accelerated flow in region III is slightly faster than 

848 in case Bioo because the pressure gradient and advection terms are 27—40 % 

849 larger in this region of Bgs- This would lead to a minimal temperature in- 

850 crease ~3 %, most notably in regions II and I (see Fig. [TOl). 

851 Outer boundary conditions for an expanded magnetosphere 

852 For our expanded configuration, we found that /poo=80 MA gave a smooth 

853 profile (least variance in FAC density). This change in /poo produced corre- 

854 sponding changes in model outputs, which are far less significant than those 

855 from our compressed case. Both magnetospheric and thermospheric angular 

856 velocities in case Cgo have slightly smaller magnitudes in the magnetodisc 
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867 region when compared to Cioo- The current profiles calculated by setting 

858 /poo=80MA produced currents strongly resembling those for case Bgs- Peak 

859 values for currents at the region III / II boundary are ~70-80 % of those 

860 in case Cioo- Thermospheric fiows for Cso differ slightly from Cioo, most 

861 significantly in the larger degree of sub-corotation in Cso- The 'hotspot' in 

862 the polar region is ~2 % hotter than in case Cioo because of faster pole- 

863 ward fiows transporting heat more efficiently. These faster fiows are due to 

864 stronger advection in case Cso producing stronger acceleration compared to 

865 Cioo- 

866 4-4- Effect of outer boundary conditions on ionospheric powers 

867 In this section, we examine figures for ionospheric power per unit area 

868 and their respective integrated power per hemisphere for cases A45, Bgs and 

869 CsO- 

870 

871 The power per unit area for case A45 in the high latitude region is shown 

872 in Fig. [TTk . Colour conventions and labels are the same as those in Figs. [Tj-l. 

873 Powers per unit area (in Fig. [TTk) are integrated over each hemisphere and 

874 are shown in Fig. [TTl l. The transition between region III and II is the most in- 

875 teresting for comparison; Aioo has a large prominent peak in magnetospheric 

876 power whereas A45 has a significantly reduced peak due to smaller values of 

877 VLm- In region II i7A/~0.5i7j (see Fig. IHk) implying that magnetospheric and 

878 atmospheric power in this region are equal. Joule heating and ion drag are 

879 thus increased in A45 compared to Aioo to meet the above requirement. These 

880 results suggest that smaller Ip^o values will generally dissipate more heat in 

881 the atmosphere but less efficiently maintain corotation in the magnetosphere. 

882 

883 For cases Bgs and Cso the power per unit area and integrated powers per 

884 hemisphere are shown in Figs. [TTb-c and e-f respectively. In comparing these 

885 two cases with Bioo and Cioo "we find only small differences in atmospheric 

886 powers (Joule heating and ion drag), predominantly at the region III / II 

887 boundary where there are two peaks with a small trough in between. The 

888 integrated atmospheric powers per hemisphere thus remain relatively uni- 

889 form with the changes in /poo specified for all baseline (Bes and Bioo) and 

890 expanded (Cso and Cioo) cases. The magnetospheric power per unit area 

891 for Bgs and Cso has significantly smaller magnitudes in region III compared 

892 to their /poo=100MA counterparts. Therefore, for these configurations of 

893 the magnetosphere, decreasing the value of /poo decreases the efficiency with 
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894 which the atmosphere can accelerate the magnetosphere towards corotation, 

895 but has no significant effect on atmospheric powers. 



5. Conclusion 
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In this study, we have expanded on the model of SA09 and described the 
effects of different solar wind dynamic pressures on the coupled ionosphere- 
magnetosphere system at Jupiter. We constructed three typical magneto- 
spheric profiles (see Tabled]), compressed, baseline (average) and expanded. 
These were then coupled to our global two-dimensional thermospheric model 
(ISmith and Aylwardl . 120081 ) and a global conductivity model of the iono- 
sphere (GG). This allowed for a comparison with results from SA09, but also 
provided a first quantitative investigation of how ionospheric, thermospheric 
and magnetospheric parameters were affected by differing solar wind condi- 
tions. 

Our results confirm many r e sults from previous studies suc h as those of 
Southwood and Kivelson ( 200l[ ): Cowley and Bunce ( 2003a ) and Cowley et al 
( 120071 ). We see an increase (resp. decrease) in thermospheric and magneto- 
spheric angular velocities for compressed (resp. expanded) magnetospheres 
relative to our baseline. The thermosphere super-corotates just equatorward 
of the middle / outer magnetosphere boundary similarly to SA09. We solve 
for Qm self-consistently in the magnetodisc in all cases using the equations 
of disc dynamics. The Qm value in the outer m agnetosphere i s a co nstant, 



dependent on disc radius i.e solar wind pressure (ICowley et al.l . l2005l ). Mag- 



netospheric angular velocities in t he polar cap, are a lso fixed at a set fraction 



10 %) of rigid corotation [Qj) (jlsbell et al.l . Il984l ). We also found that the 



coupling currents showed an increase (~20 %) in intensity when going from 
an average to a more expanded magnetosphere and a decrease (~40 %) when 
going from average to compressed. 

Our thermospheric model was used to simulate azimuthal and meridional 
neutral velocities. We see super-corotation in the azimuthal fiows equator- 
ward of the edge of the magnetodisc fiux shells. There lies a strong sub- 
corotational jet at mid to upper altitudes in the mapped ionospheric loca- 
tions of the outer magnetosphere and polar cap. The spatial size of the strong 
sub-corotation region increases with increased magnetospheric size due to the 
weaker magnetic field strength in expanded magnetospheres; thus the trans- 
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930 fer of angular momentum is less effective at maintaining corotation. Angular 

931 momentum is transferred from the thermospliere to the magnetosphere, in 

932 order to accelerate the latter towards corotation. If the thermosphere itself 

933 is significantly sub-corotating, then there is a lower 'reservoir' of available 

934 angular momentum that can be transferred. This results in a decreased 

935 plasma angular velocity in these outer regions of the magnetosphere. We 

936 see a meridional flow directed polewards at low altitudes and equatorwards 

937 at high altitudes. From the poleward edge of the magnetodisc to the centre 

938 of the polar cap, a region of accelerated poleward flow exists whose veloc- 

939 ity magnitude increases from a compressed to an expanded magnetosphere. 

940 This occurs because there is a force imbalance in this region that increases 

941 advection of momentum in expanded magnetospheres. Advection restores 

942 balance which results in the acceleration discussed above. This accelerated 

943 flow produces a 'hotspot' in the polar cap, with a maximum temperature in- 

944 crease of ~130 K from compressed to expanded magnetosphere. The size of 

945 the 'hotspot' also increases with an expanding magnetosphere. We find that 

946 the outer magnetosphere and polar cap are most strongly heated by Joule 

947 heating and ion drag. This heat is then distributed across the polar region 

948 via advection rather than viscous transport, whilst more equatorial regions 

949 are significantly cooled. This aspect of thermospheric flow is consistent with 

950 those presented in SA09. These results also suggest that accurate measure- 

951 ments of ionospheric temperature in the polar region could potentially be 

952 used to probe magnetospheric conditions. 



953 



954 We also showed that the power dissipated in the upper atmosphere (con- 

955 sisting of both Joule heating and ion drag) increases with an expanded mag- 

956 netospheric configuration. The power used to accelerate the magnetospheric 

957 plasma initially increases as we expand the magnetosphere from compressed 

958 to average configurations, but then decreases with an expansion from aver- 

959 age to expanded. This suggests that power used to accelerate the magneto- 

960 sphere has a 'local' maximum for a magnetosphere size somewhere between 

961 a compressed and expanded configuration. The total power extracted from 

962 planetary rotation is the net sum of the atmospheric and magnetospheric 

963 powers and this is positively correlated with magnetosphere size. Compar- 

964 ing our compresse d and average r nagne tospheres with the 'intermediate' and 



965 'baseline' cases in ICowley et al.l (120071 ) . we showed that the use of a two- 



966 dimensional thermosphere model results in the transfer of ~20 % more en- 

967 ergy from the thermosphere to the magnetosphere in order to accelerate the 
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968 plasma in the magnetodisc. Using our more realistic model of thermosplieric 

969 flow also produced increased dissipation of energy in the thermosphere vi a 



971 



1004 



970 Joule heating and ion drag than the cases presented in I Cowley et al.l (120071 ) 



972 We have shown that our original compressed case has some unusual cur- 

973 rent density features due to a relatively high value for the radial current at 

974 the outer disc boundary. In order to confirm this we decreased the bound- 

975 ary value of Jpoo for each case in order to produce alternative models with 

976 minimum variance in their FAC profiles. This led to the selection of Jpoo of 

977 45 MA, 68 MA and 80 MA for the compressed, average and expanded cases 

978 respectively. 

979 

980 Decreasing the radial current /poo at the boundary between the middle and 

981 outer magnetospheres resulted in all magnetosphere-ionosphere coupling cur- 

982 rents being reduced in accordance with the new value of /poo- This is expected 

983 under the assumption of current continuity. The main differences between 

984 cases with large and reduced radial currents lies mainly within the magne- 

985 todisc. For the FAC density, changes related to /poo were also significant 

986 throughout the outer magnetosphere. Thermospheric and magnetospheric 

987 angular velocities changed only slightly for the baseline and expanded case 

988 but much more substantially for our compressed case. For azimuthal flows we 

989 found that decreasing /poo also generally increased the level of sub-corotation 

990 throughout high latitudes. For meridional flows we found slight increases 

991 in localised regions of accelerated flow, most evident in the alternate com- 

992 pressed case. We also found that the polar region becomes slightly warmer 

993 with a decrease in /poo; peak temperatures for the alternative configurations 

994 increasing relative to their /poo=100MA counterparts. The total integrated 

995 powers increased with decreasing /poo for our compressed case, but decreased 

996 for our baseline and expanded cases. The integrated magnetospheric power 

997 for all cases decreased along with /poo, whilst atmospheric power increased 

998 by ~20 % for the alternate compressed case but remained almost equal for 

999 our baseline and expanded cases. Thus, it seems that decreasing the bound- 

1000 ary radial current /poo effectively decreases the 'ability' of the thermosphere 

1001 to transfer angular momentum to the magnetosphere. This behaviour, as 

1002 expected decreases the intensity of auroral emissions and produces a slightly 

1003 warmer polar region. 



1005 Our calculations suggest that main oval auroral emissions and brightness 
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1006 for an expanded magnetosphere would generally be greater than that of a 

1007 compressed one. The detailed structure of the FAC density profile in the 

1008 magnetodisc is most sensitive to the value of /poo for the compressed case. 

1009 Compressed magnetospheres in the steady state have larger field strength 

1010 than expanded ones and are more efficient at maintaining the co-rotating 

1011 magnetodisc plasma at larger distances. This leads to a smaller shear in 

1012 angular velocity between the magnetosphere and thermosphere and, conse- 

1013 quently, smaller thermospheric temperatures. As a result, auroral emission 

1014 is brightest for the most expanded magnetospheric systems. We also saw 

1015 that auroral emissions would increase at the boundary between the outer 

1016 magnetosphere and the polar region with magnetospheric compression due 

1017 to the large change in plasma angular velocity at this boundary. Better ob- 

1018 servational constraints of Qm are required to confirm this prediction. 

1019 

1020 This aspect warrants further investigation since we have not attempted 

1021 to model the change in polar cap angular velocity with solar wind dynamic 

1022 pressure. Furthermore, the caveat with these predictions is that the system 

1023 is in a steady-state (where there is no explicit time dependence of the model 

1024 outputs). We thus view this study as an initial step towards developing a 

1025 model to study the transient effects of rapid changes in the solar wind dy- 

1026 namic pressure. Results of such studies could provide further insights to the 

1027 'energy crisis' at Jupiter (SA09), and the physical origin of transient auroral 

1028 features. 

1029 

1030 Finally, the results presented in this study contribute to a larger set of 

1031 theoretical investigations which have provided useful quantitative predictions 

1032 of how the Jovian aurorae would respond to changes in solar wind dynamic 

1033 pressure. Such results are useful for interpreting auroral observations, and for 

1034 making more extensive use of such data as remote diagnostics of the physical 

1035 state of the Jovian magnetosphere. 
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Figure 1: (a) Variation of the magnetic field strength (log scale) with equatorial radial 
distance within the magnetodisc for the three configurations used. Case A is represented 
by the red solid line, whilst cases B and C are represented by the blue and green solid lines 
respectively, (b) The corresponding flux functions for the three magnetospheric cases are 
plotted against equatorial radial distance using the same colour code. 
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Figure 2: Thermospheric and magnetospheric angular velocities for cases A-C are plotted 
in the high latitude region and are represented by solid and dashed lines respectively. 
Red lines represent case A, blue case B and green case C. The black dot labelled 'lo' 
indicates the magnetically mapped position of the moon lo's orbit in the ionosphere. The 
magnetospheric regions (region III is shaded) considered in this study are labelled and 
separated by the dotted black lines. 



Table 1: Table showing the three different magnetospheric configurations used in this 
study. The radii of the magnetodisc Rmm and magnetopause Rmp are shown along 



Huddleston et al.l (|l998l) magnetopause models (J or H respectively). 



with t he values of the perturbation field. Note Rmp is calculate d as in ICowlev et aL. 
( 2007). Solar wind dyna mic pressure {Psw) is also shown for both lJov et al.l ( 2002 ) and 



Case 


A 


B 


C 


Rmm/ Rj 


45 


65 


85 


Rmp/ R-j 


75 


86 


101 


AB,/ iiT 


-1.16 


0.0 


0.19 


Pswj/ nPa 


0.121 


0.060 


0.020 


Pswh/ nPa 


0.034 


0.018 


0.008 
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Figure 3: (a) Height-integrated Pedersen conductivities for cases A-C plotted versus 
latitude. Cases A-C are represented by red, blue and green solid lines respectively. The 
magnetically mapped location of lo in the ionosphere is labelled and marked by the black 
dot. Magnetospheric regions (region III is shaded) are labelled and separated by dotted 
black lines, (b) 'Slippage' parameter K plotted versus latitude for cases A-C. The colour 
code for cases A-C remains the same as (a). 
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Figure 4: (a) Azimuthally-integrated Pedersen current shown as a function of latitude for 
cases A-C. Case A is represented by the sohd red hue, case B by the blue line and case 
C by the green line. The magnetospheric regions (region III is shaded) are also marked 
and separated by the dotted black lines, (b) Azimuthally-integrated radial current plotted 
against equatorial radial distance from Jupiter for cases A-C. The colour code is the same 
as in (a). 
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Figure 5: FAC densities in the high latitude region for cases A-C. Red sohd hues represent 
FACs for case A whilst blue and green solid lines represent FACs for cases B and C 
respectively. The magnetospheric regions (region III is shaded) are labelled and separated 
by black dotted lines. 
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Figure 6: Force balance diagrams for cases A (left column) and C (right column) at 
ionospheric co- latitude of 75 °. (a)-(b) show meridional and zonal force balance in the low 
altitude region whilst (c)-(d) show meridional and zonal force balance in the high altitude 
region. Ion drag forces are represented by blue lines, fictitious (Coriolis) forces by green 
lines, pressure gradient by red and advection by the cyan line. The velocity vector is also 
plotted and is represented by the magenta lines. Note that the magnitude of velocity 
components have been divided by a factor of 1x10'* to fit the plotted scale and that in 
(c)-(d) the components of ion drag and advection have been multiplied by a factor of 10 
to increase visibility. 
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Figure 7: Figures a-c show the variation of thcrmospheric azimuthal velocity (colour scale) in the 
corotating reference frame for cases A-C respectively. Positive values dark red represent super-corotation, 
whilst negative values (light red to blue) represent sub-corotation. The direction of meridional flow is 
indicated by the arrows and the white line represents the locus of rigid corotation. Magnetospheric regions 
(region III is shaded) are labelled and separated by black dotted lines. Figures d-f show the meridional 
velocity in the thermosphere for cases A-C. The coleur scale indicates the speed of flows. Other labels 
and lines are as for (a)-(c). Figures g-i show thcrmospheric temperature distributions. Magenta contours 
enclose regions where Joule heating exceeds 20Wkg~^, solid grey contours enclose regions where ion drag 
increases the kinetic energy at rates exceeding 20Wkg~^ and dashed grey contours enclose regions where 
ion drag decreases the kinetic energy at rates greater than 20 Wkg~^. Figures j-1 show how the power per 
unit area varies for cases A-C. The blue line represents total power which is the sum of magnetospheric 
power (red line) and atmospheric power (green line); atmospheric power is the sum of both Joule heating 
(black solid line) and ion drag (cyan solid line). Other labels are as for (a)-(c). 



Integrated power per hemisptiere for cases A, B and C 




Case 

Figure 8: Integrated ionospheric powers per hemisphere for cases A-C are represented 
in this figure. Ion drag is represented by blue bars, Joule heating by green bars and 
magnetospheric power by red bars. The white dashed line shows the division in powers 
between closed and open field line regions. Powers in the closed field regions lie below 
the dashed white line whilst powers in the open field regions lie above it. Total power 
dissipated for each mechanism (in TW) is printed on its respective colour bar. 
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Figure 9: Figures a-f show thermospheric and magnetospheric angular velocities, height- 
integrated Pedersen conductivity, 'slippage' parameter, azimuthally-integrated radial cur- 
rent, azimuthally-integrated Pedersen current and FAC density respectively for case A 
with /poo=45MA (A45) represented by red lines and case A with /poo=100MA (Aioo) 
represented by blue lines. Note that case Aioo is the same as case A in section [3] For (a) 
the solid lines represent the thermospheric angular velocity and the dashed lines represent 
the magnetospheric angular velocity. Magnetospheric regions (region III is shaded) are la- 
belled and separated by dotted black lines. Figures g-i show this thermospheric azimuthal 
velocity, meridional velocity and temperature distributions in the high latitude region for 
case A45. Arrows, contours and colour bars are the same as in Fig. [71 
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Figure 10: Figures a-f show thermospheric and magnetospheric angular velocities, height- 
integrated Pedersen conductivity, 'slippage' parameter, azimuthally-integrated radial cur- 
rent, azimuthally-integrated Pedersen current and FAC density respectively for case B 
with /poo=68MA (Bgg) represented by red lines and case B with /poo=100MA (Bioo) 
represented by blue lines. Note that case Bioo is the same as case B in section |31 For (a) 
the solid lines represent the thermospheric angular velocity and the dashed lines represent 
the magnetospheric angular velocity. Magnetospheric regions (region III is shaded) are la- 
belled and separated by dotted black lines. Figures g-i show this thermospheric azimuthal 
velocity, meridional velocity and temperature distribution in the high latitude region for 
case Bgs. Arrows, contours and colour bars are the same as in Fig. [71 
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Figure 11: Figures a-c show ionospheric powers per unit area in the high latitude region 
for cases A45, Bgs and Cgo respectively. Total power per unit area is represented by the 
blue line, magnetospheric power by the red line, atmospheric power by the green line. 
Joule heating by the black line and ion drag power by the cyan line. Magnetospheric 
regions are labelled and separated by dotted black lines. Figures d-f show integrated 
powers per hemisphere for cases A45, Bgs and Cso respectively. Ion drag. Joule heating 
and magnetospheric powers are indicated by blue, green and red bars. Powers in the closed 
field line regions lie below the white dashed line whilst powers in the open field regions lie 
above it. Total power dissipated for each mechanism (in TW) is printed on its respective 
colour bar. 
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